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ABSTRACT. We give an application of the Crandall-Rabinowitz theorem on local bifurcation to a system of 
nonlinear parabolic equations with nonlocal reaction and cross-diffusion terms as well as nonlocal initial con- 
ditions. The system arises as steady-state equations of two interacting age- structured populations. 



1. Introduction 

In this note we consider coexistence solutions to age-structured population dynamics with diffusion, the 
main feature of the present work being the inclusion of nonlocal cross-diffusion terms. More precisely, we 
shall establish positive nontrivial solutions u = u(a, x), v = v(a, x) to the system 

d a u — div x (d\iy)V x u + uX7 x d,2(V)) = — au 2 — fj,i(V)u , a e (0, a m ) , x e f2 , (LI) 

8 a v -div x (d 3 (U)V x v + v\7 x d 4 (U)) = -Pv 2 - fi 2 (U)v , ae(0,a m ), x e , (1.2) 

for a e (0, a m ), and x e O, subject to the nonlocal initial conditions 

u(0,x) = r]U, ie!l, (1.3) 

v(0,x) = £V, (1.4) 

and Dirichlet boundary conditions 

u(a, x) = , ae(0,a m ), x e dfi , (1.5) 

v(a, x) = , a e (0, a m ) , x e dfl , (L6) 
where we agree here and in the following upon the notation 



U 



ui(a)u(a, ')da , U := b(a)u(a, •) da (L7) 
Jo Jo 

for the function u defined on J := [0, a m ] and analogously for the function v. Equations dl.U - dl.7t arise 

naturally as steady-state (i.e. time-independent) equations of two age-structured populations with densities 

u and v, respectively, and maximal age a m > living in a (bounded and smooth) domain fl <z K™, where a 

is the age and x is the space variable. The integrals with respect to age in (ll.lt and dl.2t are (weighted) local 

total populations with a given nonnegative weight function u>. The divergence terms in dl.lt . dl.2t describe 

spatial movement with nonlocal coefficients dj . They reflect intrinsic dispersion as well as an increase of 

dispersive forces by repulsive or attractive interferences with an increase of the other population. We refer 

to lfl2l for a derivation of such kind of models (without age-structure). The right hand sides of dl.lt and 

dl.2t take into account intra- and inter-specific interactions of the two populations with constants a, j3 > 

and functions depending nonlocally on the population densities. Creation of new individuals is described 

by (II. 3t . (II. 4t with birth profile b and parameters r], £ measuring the intensity of the fertility. We reference 

to ifTsl and lfl5l for further information on the modeling assumptions. To avoid unnecessary notational 

complications, the equations above are stated as a simplified version of more elaborate models, and we 

remark that the subsequent analysis would not change in any way if one would allow for e.g. different 

weight functions, different birth rates, or different maximal ages for the two populations. 
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In this note we give an extension of previous results fl31 [TBI [17 1 being described in more detail in the 
next section. Steady-states for a single age-structured populations were investigated e.g. in [8|. We shall 
also point out that steady-state solutions for two interacting populations when age-structure is neglected, 
i.e. variants of the elliptic counterparts of (ll.l| )-( fL2] i, have attracted considerable interest in the past, see 
for example O [5] |7] [9] [TO] [H] [T31 and the references therein. 

2. Notation and Main Result 

The main features of the equations under consideration are the nonlocalities appearing in the diffusion 
and reaction terms as well as in the initial conditions. Similar equations with local reaction terms (i.e. 
l-ii(v) = a^v, H2(u) = with Qf2,/32 > 0) have been investigated in ffT31 [TBI with linear diffusion 

(i.e. d\ = d,3 = 1 and d,2 = g?4 = 0) and in ifTTl with a local cross-diffusion term (i.e. di(v) = 1 + v, 
d2(v) = v, c?3 = 1, and = 0). In these papers global bifurcation results have been derived with respect 
to the parameters rj and £. The aim of this note is to show that (local) bifurcation results can be obtained for 
equations (11.1 b - (fTT7t including nonlinear nonlocal diffusion terms. More precisely, we shall provide values 
for the parameters r\ and £ for which dl.U -( [L7b have coexistence solutions, i.e. smooth solutions (u, v) with 
both components nontrivial and positive. Establishing positive steady-state solutions is a first step toward 
an understanding of (time-dependent) two population dynamics. To this end, we assume throughout the 
paper that 

di, da^^^eC^R), d 2 , d 4 e C 2 (R) , (2.1) 

satisfy 

d 2 (0) =/ii(0) = 0, (2.2) 

respectively, 

dj(z)^S, z^O, j = l,3, (2.3) 
for some 8 > 0. For an easier reference in the future we suppose that 

di(0) = l. (2.4) 
For the weight and the birth functions we assume 

w,6eL+(J), 6(a) > for a near a m (2.5) 

together with the normalization 

ra m 

b(a)e- Xia da = 1 , (2.6) 

Jo 

where Ai > denotes the principal eigenvalue of — Aa; on fl subject to Dirichlet boundary conditions. For 
technical reasons we introduce h q := L q (J, L q (Q)) and the solution space 

W q :=L q (J,Wl D )nW^(J,L q ) 

with q e (n + 2, oo) fixed, where L q := L q {Q), and W* D := W£ D (£l) refer to Sobolev-Slobodeckii spaces 
including Dirichlet boundary conditions if meaningful, i.e. if k > 1/q. We let W+ denote the positive cone 
of W q and set W+ := W+\{0}. RecaU the embedding 

W^C 1 - 1 ^-*^,^],^), 0<1?<l-l/ ? , (2.7) 

which holds for § = 1 — 1/q due to 12 III. Thm.4. 10.2] and otherwise by the interpolation inequality ||2] 
I. Thm.2. 11.1]. In particular, the trace 70 u := u(0) defines a bounded linear operator 70 6 £(W 9 , W 2 ~^ q ). 
Also recall (e.g. from [2|) that A e L Vj (J, £(W q 2 D , L q )) is said to have maximal L q -regularity provided 
that the operator 

(d a +A, l0 )e£(W q ,h q xW q 2 - D 2/q ) 
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has a bounded inverse, where (A(f))(a) := A(a)4>(a) for a e J and </> e W 9 . We also note the compact 
embeddings 

KD^rf- 2/q ^C\n) (2.8) 

as q > n + 2, and that the interior of the positive cone W 2 ^^ q ' + of W 2 ^^ 9 , denoted by mt(W 2 ^ q ' + ), is 
nonempty. Finally, recall that an operator Z e C{W 2 ~^ q ) := C(W 2 ~^ q , W 2 ~^^ q ) is strongly positive if 
e mt(W^ 2 *) for e < D 2/9 ' + \{0}. 

On taking v = 0, our assumptions imply that we get from dl.lt -( [TT7l i a reduced problem with a nonlocal 
initial condition of the form 

d a u - A x u = -au 2 , u(0) = rfU , (2.9) 

subject to Dirichlet boundary conditions, which has been studied in a previous work lfl5l Thm.2.1]. In fact, 
in view of the normalization (12.61 . this problem admits no solution in W+ if r\ ^ 1, while for each r/ > 1 
there is a unique solution u v e W+ depending smoothly on n and ||urj||w, — > oo as r\ — > oo. Parabolic 
regularity theory implies that «„ is smooth with respect to both a e J and x e il. 

In the following we shall consider the situation that r\ > 1 is given. We then regard £ as a bifurcation 
parameter and write (£, u, v) for solutions to dl.li - dT?7Y The considerations above guarantee the existence 
of a semi-trivial branch of solutions 

To :={(£,«,,<>); £>0}. 

On applying the famous Crandall-Rabinowitz theorem on local bifurcation, we derive the existence of a 
local branch of coexistence solutions bifurcating from the semi-trivial branch To. More precisely, we have: 

Theorem 2.1. Given rj > 1, there exists £o := £o(^) > smc/; f/zaf a /oca/ curve T in M + x W+ x W+ 
of (smooth) coexistence solutions (£,u,v) to ( 1 1 . lb - dTTTb emanates from (£o,w^,0) 6 T> bifurcation 
point is unique, i.e. there is no other bifurcation point on To to coexistence solutions. 

As equations d 1 - 1 b - dT/Tb are symmetric in u and v, one may interchange the role of £ and n, of course. We 
remark that establishing a bifurcation from the other semi-trivial branch {(£, 0, ug); £ > 1} when regarding 
£ as bifurcation parameter (if /X2(0) = and ^(0) = 0) does not seem to be obvious for the present 
situation. However, for the case of linear diffusion, i.e. if di = da = 1 and d 2 = d± = 0, one can show 
global bifurcation results for jl.lMl.7i along the lines of lfl5l[T6l for nonlocal reaction terms as well. 



3. Proof of Theorem [271 



The remainder is dedicated to the proof of Theorem l2.ll which is a consequence of the bifurcation result 
of Crandall & Rabinowitz [4, Thm.1.7]. Let r\ > 1 be fixed. We write solutions to dl.li - dl.7i in the form 
(£, u, v) = (£,,u v — w, v), which we then obtain as the zeros (£, w, v) of the function 



F 



h q xh q x W 2 -*l q x W 2 ~ D 2/q 



where 



( 3 a w - div x (d 1 (V)V x w + wV x d 2 (V)) - div x ((l - dxiV)) d2(V)) \ 



d a v - div x (d 3 (U v 



— — w) + loai^w — 
W)V x v + vV x d A (U v - W)) + Pv 2 + fi 2 (U v 
w(0) -rjW 
v(0) - £V 



W)v 
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for (£, w, v) e K x W q x W 9 . Owing to d272b and d274b , the Frechet derivatives at (£, to, u) = (£, 0, 0) are 
given by 



/ <9 a - A x (j) + div x (di (0)* V 4(0)V^) -/4(0)*M r) +2au v 4> \ 

d a ip - div x (d 3 (U v )V x ip + ipV x d4,(U v )) + ^2(^)1/' 
0(0) - 77$ 

^(o) - 



with dashes denoting derivatives, and 

^,(»,t,)(e,o,o)[^^] 



/ \ 





V -* / 

where we use here and in the following notation dl.7| > for (0, tp) e W q x W q . 



(i) We shall show that F/ WiV ^(£, 0, 0) is a Fredholm operator of index zero for £ > 0. To this end we 
introduce operators A\ £ C 1 / 2 ( J, C(W^ D ,L q )) and A3 e Z^W^ , L g ) as 

-Ai(a)0 := — + 2au v (a)(f> , 

:= -div x (<i3(^ I7 ) V ^+^ V » d 4(L\)) + ^(Ur,)^ , 

for a e J and 4>,ip e W qD , the regularity being implied by (I2.71 i as u,, e W q . It follows from (12.31 i and J2] 
III.Sect.4] (in particular, see III.Ex.4.7.3d), III.Thm.4.8.7, and III.Thm.4. 10.10 of 0) that both A x and A 3 
have maximal L g -regularity. We also introduce A2 e C(W q , L g ) by 

(A a ^)(a) := div a: (di(0)*V a! u,(a) + ii„(a)d' 2 (0)V x *) -/ii(0)^(a), a e J , ^eW,, 
and set, for (0, -0) e W 9 x W,, 

' A l <p + A 2 1p 



A(0,V) 



using the convention (Ai4>)(a) = Ai(a)0(a) and (^^(a) = A 3 ij}{a) for a e J. On letting 



7o (0» 



0(0) 
<0(O) 



,V>) 6 Ei 



and 



Eo:=L,xL,, B f := W*^ 2 /<? x W^ D 2/q 



it is readily seen from the triangular structure of A that the operator (d a + A, 70) 6 £{Ei,Eo x£ 5 ) has a 
bounded inverse, say, 

T:= (5 a +A,7o) _1 e£(E x£ c ,Ei). 



Introducing £[£] e £(Ei,£i) for £1 := W 2 ^ x by 



^](0,^) 



77$ 



we obtain from a straightforward modification of [ 14, Lem.2.1] (with Eq := L q x L q therein) that 

:= ^(»,,)K,0,0) = (d a + A, 7 o 6 £(Ei,E x 

is indeed a Fredholm operator of index zero. In fact, defining 

Q[€\w:=e[S](T(0,w)) , W e£ s , 
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and observing that Q[£] e £(E^,Ei) is a compact operator on E q into itself due to the compact embedding 
Ei «-> E$ implied by ( 12.8b . we have 

ker(LK]) = {T(0,w) ; w 6 ker(l - Q[£])} , 

rg(L[£]) = {(/, h) e E x E, ; ft + £[£](T(/, 0)) e rg(l - Q[£])} , 

both spaces being closed with 

dim(ker(L[£])) = codim(rg(L^])) = dim(ker(l - Q[f])) < oo . (3.1) 

(ii) We now choose £o > such that ker(L[£o]) is one-dimensional. First observe that, owing to the 
parabolic maximum principle (e.g. see E Cor.13.6]), the semigroup {e~ aA3 ; a ^ 0} on L q generated by 
—A% is such that e~ aAs e C(W^~^^ q ) is strongly positive for a > 0. Hence, (12.5b and standard regularity 
effects of analytic semigroups [2| ensure that 

H := ^ 6(a) e~ aA > da e C{W^ a/q , W^) 
Jo 

defines a strongly positive and compact operator on W^~^ q into itself by ( 12.81 ). Its spectral radius r(H) > 

is thus a simple eigenvalue of H with a corresponding eigenfunction e iat(W^^^' + ) due to the Krein- 
Rutman theorem [ 1 , Thm.3.2]. We then put 

and obtain ker(l — t;oH) = spanj^o}- Let us also observe that A\ e C 1 / 2 ( J, C(W^ D , L q )) generates a 
parabolic evolution operator Ui(a, a) on L q by 12 Cor.4.4.2]. The same arguments as above ensure that 

Gi := f m 6(a)[/ 1 (a,0)da6/:(^ 2/9 ,^ 2 I) ) , 
Jo 

that Gi is compact on ^ q into itself by (12.81 l. and that Gi is strongly positive by ( 12.5b since this is 
true for each Ui(a, 0), a 6 (0, a m ], see Cor.13.6]. Next we claim that 1 — r\G\ is invertible. Indeed, 
owing to the positivity of u v , the evolution operator U2(a, a) on L q generated by At — au v = —A x + au v 
(subject to Dirichlet boundary conditions) is such that [^(d, 0) — U±(a, 0) is strongly positive on W q 2 ~^ q 
for a e (0, a m ] by the parabolic maximum principle, whence G2 — Gi is strongly positive due to ( 12.5b . 
where 

G 2 := f b(a)U 2 (a,0)da . 
Jo 

Therefore, r(Ga) > r(Gi) by the Krein-Rutman theorem JT] Thm.3.2]. Recall then that ( 12.9b implies 
u v (a) = U2(a,0)u v (Q) with u^(0) = rjG2U rj (Q) and so r^G-z) = 1 again by the Krein-Rutman the- 
orem since u v (0) is strictly positive. Consequently, 1 — r\G\ is invertible. Consider now the equation 

£[£o](0, VO = for (0, 0) e Ei, that is, 

d Q + A 3 V> = , a e (0, a m ] , ^(0) = £0* , (3.2) 

and 

d a <l> + A^aty = -(A 2 ip)(a) , ae(0,a m ], 0(0) = r?$ . (3.3) 

From (O we deduce -0(a) = e~ aA3 ip(0) with -0(0) = ^ o Hi/}(0), whence -0(0) = C*o for some (el 
since 1 is a simple eigenvalue of £qH with eigenvector ^q, and so ip = (ip* with 0*( a ) : = e ~ a ^o- 
From the first part of ( |3.3b it then follows that 

0(a) = Cfi (a, 0)0(0) -C f r7i(o,o-)(i4 a ^)(CT)d<T, aeJ, 
Jo 
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which, when plugged into the initial condition (f)(0) = 77$, yields (f>(0) = (&o, where 

$0 := -il(l - ndy 1 \ b(a) \ U 1 (a ) a)(A 2 ifj*)(a)dada. 
Jo Jo 

Thus, setting 

r a 

<^(a) := J7i(a,0)$ - Ui(a, a) (A 2 ip^)(a) da , aeJ, 
Jo 

we conclude ((f>, ip) = C(4>*ii } *) an d then 

ker(L[ft]) = span{(0„^,)} . (3.4) 

(iii) Next, we check the transversality condition of |4|, that is, we show that 

i ? C,( U)lO )(&,O,O)[^ ) ^]0rg(LKo]) ■ (3.5) 
Supposing otherwise, there is ip e W q such that 

d a ip + A 3 ip = 0, ae(0,a m ], ip(0) - ft* = , 
and we easily derive (1 — fti/)-0(O) = — Vt*. Since 

*o e int(W^ 2/9, + ) n ker(l - ftff) , 
we may choose n > sufficiently large such that 

P := «tf - i/>(0) e mt(W*- 2/q > + ) , (1 - ftF)P = * t . 

Because r(fti/) = 1 and "J* e W^~^ q ' + \{0}, this contradicts the fact that this last equations has no 
positive solution P according to 0] Thm.3.2], and we conclude ( 13. 5t . 

(iv) Gathering ( 13. It . ( 13.4b . and (13.51 ) we are now in a position to apply [4 Thml.7] and deduce that 
the nontrivial zeros of F close to the bifurcation point (ft, 0, 0) lie on a continuous curve. More precisely 
and applied to dl.U - dl.7b . we obtain some £0 > and continuous functions £ : (— £o,£o) — * K an d 
Cj : (— £q,£o) — * W 9 with £(0) = ft and Cj'(O) = 0, such that the nontrivial solutions (£,u,v) to dl.lt - 
dl.7t lie on the curve 



{(£(£),«»? -£Ci(£),#* + eft(e)); -e < e < £0} <= 



'9 > 



which bifurcates from (ft, u v , 0). Since w^(0) and (0) = 'So belong to int^WJ 2 ^, ' + ), it follows from 
d2.7l i that, for g > sufficiently small, u £ := it r) — £0* — eft (e) and w e := ei/i* + eft (e) both have nontrivial 
initial values 

u"(0) = u n (0) - e$ - £7oCi(e) e W*j? /q ' + 

and 

t; e (0) = e* + £7oft(e) e <,L 2/9 ' + , 
respectively. Therefore, making £0 smaller, if necessary, we derive from dl.ll ). (1 1 -2b , and the parabolic 
maximum principle that 

T := {(fte),w,, - e</>* - eftXs),^* + eft^)) , < e < e } 

is a curve of solutions to dl.U - dl.7t in M + x W+ x W+ bifurcating from (ft, u v , 0) e To. 

(v) Finally, if (£, 1^,0) e To is any bifurcation point to positive coexistence solutions of dl.U - dl.7t . 
there is a sequence (ft , Uj , Vj ) in R + x W+ x W+ of solutions to dl.lt -( fk7] i converging to this point. Since 
solutions to dl.lt - dl.7t are smooth by standard parabolic regularity theory, it is not difficult to see that a 
subsequence of ipj := Vj/\\vj ||w, converges to some ip e W+ satisfying 

d a ^ + A 3 i) = 0, 6(0, 0,,,], 



A NONLOCAL NONLINEAR REACTION-DIFFUSION MODEL 



7 



what readily implies £ = £o since ^(0) = £Htp(0), i.e. r(£H) = 1 by the Krein-Rutman theorem. Thus 
(Co j u n , 0) is the only bifurcation point on To- This proves Theorem l2.ll 
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